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We explore the possibility that quantum cosmology considerations could provide 
a selection principle in the landscape of string vacua. We propose that the universe 
emerged from the string era in a thermally excited state and determine, within a 
mini-superspace model, the probability of tunneling to different points on the land- 
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scape. We find that the potential energy of the tunneling end point from which the 
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universe emerges and begins its classical evolution is determined by the primordial 
temperature. By taking into account some generic properties of the moduli potential 
we then argue that the tunneling to the tail of the moduli potentials is disfavored, 
that the most likely emergence point is near an extremum, and that this extremum 
is not likely to be in the outer region of moduli space where the compact volume is 
very large and the string coupling very weak. As a concrete example we discuss the 
application of our arguments to the KKLT model of moduli stabilization. 
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I. INTRODUCTION 

The existence of a multiverse of solutions to string theory makes it eminently desirable to 
find a dynamical principle that would select a universe, or a subset of the multiverse that has 
properties similar to ours. Such a principle might be provided in the framework of quantum 

cosmology HQS 

. The wave function of the universe yields a probability distribution for 
the dynamical parameters of the multiverse. Depending on the boundary condition chosen 
for the wave function, the probability distribution is sharply peaked at a zero value of the 
cosmological constant (CC) j^J or is peaked at a large value of the CC jsl 0| giving rise to 
an inflationary universe. A related and complementary approach is based on the statistics 
of solutions , one of whose outcomes is that the number density of solutions is 

uniform as a function of the value of the CC. 

In this paper we propose that the universe emerged from the string era in a thermally 
excited state above the Hartle-Hawking vacuum. We show that imposing this boundary 
condition on the wave function of the universe leads, within a mini-superspace model, to 
interesting restrictions on the allowed dynamics. High temperature effects were first intro- 
duced in a related context, to the best of our knowledge, by Vilenkin ,6] who suggested 
including them in the context of the so-called tunneling boundary condition, but did not 
work out the consequences. We will show in detail the relevance of these effects both to the 
Hartle-Hawking (HH) wave function and to the tunneling wave function. 

Our proposal is motivated by an interesting observation that was made recently by 
Sarangi and Tye [ill (see also and [3]). They considered the tunneling amplitude 
as a function of the CC A, and found that if rather than evaluating it from the HH wave 
function 2] that describes tunneling from nothing, 

~ e 3n/GA } (1) 

a modified wave function is used, then the tunneling amplitude 

^ ^ e 3n/G\-(a/2)(3n/G\) 2 ^ /g\ 

has a critical point (a maximum for the HH choice). The additional parameter a is a 
constant that depends on the string mass scale M s , the Planck scale M p = G^ 1 ^ 2 and the 
number of fluctuating degrees of freedom n do f. The critical CC that Sarangi and Tye found 
is A c = An do} ($) Ml 
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The original argument by Hawking was that since the wave function ((TJ is peaked at 
A — ► + it explains why the observed value of the CC today is zero. Of course at that 
time there was just an observational upper bound on the CC, and it was assumed that the 
upper bound implies that the true value is zero. The problem with this argument is that 
it predicts a large but empty universe. Sarangi and Tye argue that their modified wave 
function (J2J) predicts inflation and hence makes the choice of the HH boundary condition 
more interesting. 

We use a mini-superspace model for our discussion. In this model only the time- dependent 
scale factor of the universe, and some time-dependent but homogeneous fields are considered 
and all other degrees of freedom are ignored. We regard the mini-superspace calculation as 
a toy model which by itself is a well defined quantum mechanical problem that can be solved 
self-consistently and hopefully incorporates some features of the complete quantum gravity 
calculation. At the current stage of development an attempt to do a more complete quantum 
gravity analysis seems to lead to the notorious problems of Euclidean quantum gravity. Thus 
our calculation should be regarded as an indication of the true solution of the full dynamical 
problem. For self-consistency we apply the quantum tunneling picture in four dimensional 
mini-superspace only for large scale factors and energies and temperatures that are small 
compared to the string scale and the Kaluza-Klein scale. The existence of the primordial 
temperature allows us to do this at the price of not being able to say anything about the 
quantum origin at zero scale factor. We regard the inclusion of primordial radiation as a 
parametrization of our ignorance of the physics at the quantum origin. 

We start our discussion in the simpler case for which the moduli, the dilaton S, and the 
volume modulus T are fixed by stringy effects. In this case the thermal boundary conditions 
lead to results that are similar to those of Sarangi and Tye. Then we take into account the 
dynamics of the moduli by assuming they are stabilized as in the recent works on moduli 
stabilization by fluxes and non-perturbative effects 3, 3|- Moduli stabilization by fluxes 
was considered previously in [ly, [jjl 18, Ell Q We show that in order to have a finite 
potential barrier and therefore a large tunneling amplitude, the end point of the tunneling 
should be to a region of the moduli potential that supports accelerated expansion. Taking 
into account additional features of the moduli potentials we further show that the tunneling 
takes place to a maximum or a saddle point that is within a bounded region in moduli space. 
In particular, we find that the tunneling to the tail of the moduli potentials, where they 
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runaway to decompactification for example, is disfavored. The reason for this is that these 
regions do not lead to sustained accelerated expansion. 

The organization of the paper is as follows. After a brief review of the wave function 
calculus in sect. |H]we discuss in sect. IIHI thermal effects and their significance to the tun- 
neling amplitude. In sect. IIVI we discuss the dependence of the tunneling amplitude on the 
parameters and conclude that we must add the dependence on moduli. We follow our own 
conclusion and in sect. 0and sect. IVII where we extend the analysis to a case with dynam- 
ical moduli. In sec. I VI II we discuss in detail the application of our results to the model of 
Kachru, Kallosh, Linde and Trivedi (KKLT). Section fVIIII contains our conclusions, and in 
the appendix we compare our work to that of Sarangi and Tye. 



II. REVIEW OF THE WAVE FUNCTION CALCULUS 

Let us first briefly review the calculation of the wave function of the universe in the 
mini-superspace context, Consider a theory with a dynamical CC A whose action is 

S= 1^gJ rf V^-2A). (3) 

Here all other dynamical fields are ignored and it is assumed that their effects are incor- 
porated into a single parameter A that is essentially an integration constant. The problem 
is further simplified by considering only homogeneous and isotropic metrics that describe a 
closed universe. For such metrics the line element is 

ds 2 = a 2 (-dt 2 + a 2 (t)dnl) , (4) 

a 2 = 2G/3tt. For this class of models the action simplifies considerably, 

S = ^ J dt(-aa 2 + a- Aa 3 ). (5) 

All quantities in the action (jSJ) are dimensionless, in particular, the dimensionless CC A = 
(7 2 A/3 = 2GA/97T. The canonical momentum conjugate to a is U a = —aa and the classical 
Hamiltonian constraint is 

H=~ (K + U(a)) « 0. (6) 

The last relation means that the Hamiltonian vanishes on the space of classical solutions. 
The potential U(a) = a 2 — Aa 4 is positive for < a < 1/vA and negative for a > 1/V\. 
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The quantum equation derived from the Hamiltonian constraint is the Wheeler-deWitt 
(WDW) equation which replaces the Schroedinger equation in quantum gravity. The mo- 
mentum operator fl a = —id /da replaces the classical momentum, 

HV[a] = (fl 2 a + 17(a)) *(o) = 0. (7) 

The WDW equation can be solved in the WKB approximation, 

*±[a] ~ e ±iSc = e ±if ° nada 

The last equality is valid only for < a < 1 / vA. In eq. (JBJ) we have ignored the prefactors 
because they will not be important for the rest of the discussion. The classical action S c is 
evaluated on a classical solution so that IT a = 8S C / da = \J—U(a). For < a < 1/ y/X, n a is 
of course imaginary and the wave function is real. Hartle and Hawking have imposed the 
no-boundary boundary condition which in this context means that the positive sign should 
be taken in eq. (JHJ). Then the under the barrier wave function near a = 1/ y/X is given by the 
growing exponential, and the probability of the universe emerging after tunneling through 
the barrier is given by 

Phh[X] ~ e 2 / 3A . (9) 

Hawking 4] has argued that this distribution explains the observed vanishing of the CC. 
At that time it was assumed that the observations indicated that it was in fact zero. Now 
there is strong evidence that the expansion of the universe is accelerated and therefore that 
the CC, or a similar form of dark energy is small and positive. However, as many have 
pointed out, the probability distribution (jHJ) also predicts an empty universe - there would 
have been no primordial inflationary stage. 

The result (JOJ) depends very sensitively on the nature of the state and in particular on 
the boundary conditions. If one uses the Euclidean path integral to define the state, the 
no-boundary proposal gives the ground state wave function. The resulting distribution is 
rather surprising. It implies that the tunneling amplitude increases exponentially when 
the barrier becomes larger. This is certainly not what happens in laboratory tunneling 
experiments which show that the tunneling amplitude decreases exponentially with the size 
of the barrier. 
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Linde jjj and Vilenkin [jj have proposed different boundary conditions that yield results 
that are similar to the probability distributions in standard tunneling processes. Vilenkin 
has proposed that only an outgoing wave should exist in the Lorentzian region for a > 
1/VA. This means that the universe has only an expanding component whereas the HH 
wave function is a superposition of expanding and contracting universes. The probability 
distribution that results from Linde's and Vilenkin's proposals is 

Pl,v[\] ~ e" 2 / 3A . (10) 

This clearly favors a large CC. In fact it favors a situation for which the barrier is as small 
as possible, and even no barrier at all! Since the semi-classical theory is valid only up to 
some cutoff scale, which we will choose to be the string scale M s , this seems to im ply that 
the universe is created in a state with string scale CC. Linde has proposed recently [22] that 
a flat compact universe may perhaps be more likely, since in this case there is no barrier. 
To discuss inflation it is necessary to include at least one inflaton field. This is done for 



instance in 



23] where it is shown that the formulae (j9H0|) remain valid when one replaces 
the CC with the scalar potential V(<f)), A — > 4ir 2 a 2 V((j)), provided that it is slowly varying, 
i.e. \V~ 1 dV/d<f)\ <C 1. Obviously, these semi-classical considerations are valid only for regions 
in field space where |V(0)| <C Mf. In various possible potentials are illustrated with 
the corresponding probability distributions for the two cases fl9|l()j1 . In all cases it is clear 
that the HH distribution will favor tunneling into the lowest positive points in the potential, 
while the LV distribution will favor tunneling into the highest points of the potential allowed 
by the cut off 
Vilenkin 



23] has argued that the above considerations lead to a 'prediction' of inflation 
from the LV wave function, in contrast to what is obtained from the HH wave function. 
However the problem is that the conclusion is cutoff dominated. Besides, the natural value 
of a cutoff would be close to the Planck scale, perhaps one or two orders of magnitude below, 
and this would be too high to agree with observations. The predicted Hubble parameter 



during inflation is H ~ y Mf /M 2 . If the string scale is taken to be an order of magnitude 
below the Planck scale - which is the case for string compactifications where there are no 
anomalously large extra dimensions, then H ~ 10 16 GeV, in conflict with the WMAP result 
that H < 10 l4 GeV. To get agreement, the cutoff scale would have to be at least two orders 
of magnitude below the Planck scale. 
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We would like to make two remarks about the possible application of the wave function 
calculus to string theory. First, in the string theoretic context V(<p) is typically steep, so 
that the condition ^^dV / d(f)\ <C 1 is obeyed only in the vicinity of the critical points of 
the potential. Thus, strictly speaking the analysis can only compare the relative probability 
of tunneling to different critical points. The LV wave function would predict tunneling 
to the highest critical point whilst the HH wave function would favor the lowest positive 
one. Second, recent work has shown that in string theory the acceptable solutions that 
have potentials that can stabilize the moduli necessarily involve fluxes of RR and NS-NS 
fields. The fluxes make the parameters of the potential, and in particular the CC, discrete 
integration constants. In the string theory context one should really consider a probability 
distribution that depends on all the moduli as well as on the flux parameters. In contrast, 
the original Hawking argument was made in a model without a scalar field. The CC was 
treated as a dynamical quantity arising as an integration constant characterizing the flux of 
a four form field strength. 



III. THERMAL EFFECTS 



Our boundary condition proposal for the wave function of the universe can be stated 
as follows: The universe emerges from the string era in a thermal state above the Hartle- 
Hawking vacuum. We propose that the decay of all the string excited states has created a 
primordial thermal gas of radiation at a temperature that is somewhat below the Hagedorn 
temperature. At this high temperature there would be in addition to the massless states 
some population of a Boltzmann suppressed massive string states, and perhaps also a gas 
of branes j^. These may behave as pressureless matter, or have some other behavior. 
We will ignore such contributions for simplicity, since we do not expect this to change the 
qualitative behavior that we find. Within the context of string theory the effective field 
theory arguments that we use make sense only at the end of the string era when the energy 
densities in the universe are somewhat below the string scale. 

Rather than evaluating the no-boundary thermal partition function, we evaluate the 
Euclidean mini-superspace path integral with a modified effective potential that includes 
the temperature corrections. The leading temperature correction is very simple, a negative 
term proportional to the fourth power of the temperature is added to the potential in the 
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WDW equation. Equivalently, we can solve the WDW equation with a modified potential, 
and even though one does not expect a coherent wave function to describe a thermal state, 
the square of this WDW 'wave function' would have the interpretation of a density matrix 
that measures relative probabilities. 

Since the temperature scales as the inverse scale factor the energy density of the radiation 
will be of the form Prad — K/(a 4 a 4 ) where K is a constant. The effective potential that 
goes into the WDW equation becomes 



The potential barrier is now in the region limited by the roots of U(a)—K = 0, a_ < a < a+, 



Of course to have a barrier at all the radiation term cannot be too large 4KX < 1. This 
condition needs to be satisfied for the semi-classical theory to remain valid as argued below. 

To keep the semi-classical effective field theory approach self-consistent the highest radia- 
tion density in the region of interest should be less than string scale, i.e. p(a_) ~ K/a A at_ = 
CM 4 with C < 1. Using the value a 2 _ ~ 1/2A from eq.(JT2J) we get K ~ -^M 4 a A . Expressing 
C in terms of the number of degrees of freedom n^of in thermal equilibrium , C = ri^of / c 4 
we have our final expression for the radiation energy density, 



The constant c must satisfy c 4 > n^o/ for the consistency of these arguments. This energy 
density corresponds to an initial temperature at a = a_ of T ~ M s /c i.e. a temperature 
close to the Hagedorn temperature. The same estimate can also be obtained by requiring 
that the initial entropy is close to saturating the entropy bound as discussed in the appendix. 
The condition for the existence of a barrier is thus equivalent in this context to the condition 
that the cosmological constant be smaller than the string scale. 

Including the contribution of the thermal radiation energy density into the effective action 
it becomes, 



U{a) — K = a 2 — Xa 4 — K. 



(11) 




(12) 




(14) 



where 



v = n dof 



1 M 4 
9vr 2 c 4 M 4 ' 



(15) 



9 



The corresponding classical Hamiltonian that replaces the one in eq.® is 

ff=~(n£ + tf(a)-£)«0. (16) 
The boundaries of the barrier can be rewritten as 




4 = 777 U±\ 1-T- ■ ( 17 ) 



It is important that for the range A < 1, < 4i//A < 1 the scale factor is large under 
the barrier a± ^> 1, so the semi-classical mini-superspace calculation is self-consistent. Also 
since one expects ridof ~ 10 2 — 10 3 , the above restriction on c means that it is about 10 
but need not be much bigger. In fact it is reasonable to expect that the initial temperature 
at the beginning of the classical evolution T ~ M s /c is close to but not quite at the string 
scale. The initial volume of the universe as it emerges from under the potential barrier into 
the classically allowed region is a 3 a 3 f ~ (j 3 /A 3 / 2 . The horizon volume on the other hand 
is approximately H~ 3 ~ a 3 /A 3 / 2 . Thus the universe starts its classical evolution having a 
volume which is approximately one horizon volume. 

Now, in addition to the Lorentzian region to the right of a+, there is also a Lorentzian 
region to the left of a_ . This is depicted in Fig. [3 The solution of the WDW equation in 
the WKB approximation is obtained by matching the solutions on the boundaries of the 
three different regions, 

*j = A+e +m ' + A-e'** 1 , (18) 
\S? n = B+e + *" + B^e-^", (19) 
m m = C+e +i011 ' + CLe-^'w. (20) 

The exponents are given by 

$ 7 (a) = VXJ v 7 ^- -a 2 )(a\ - a 2 ), < a < a_ (21) 
o 

a 

$//(a) =V\ J a/(o 2 - a 2 _)(a 2 + - a 2 ), a_ < a < a+ (22) 
$ m (a) =V\j ^(a 2 -a 2 )(a 2 - a 2 + ) a+ < a. (23) 



In practice the matching has to be done with care since the boundaries of the three regions 
are turning points where E = U, and therefore the WKB approximation breaks down there. 
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The new Lorentzian region a < a_ can clarify and resolve the debate about which linear 
combinations to take inside the forbidden region, and which boundary conditions to choose. 
If one puts boundary conditions in this region that correspond to "initial conditions" and 
not to "final conditions" about the state of the universe when it starts its classical evolution 
after tunneling then any generic choice will effectively be equivalent to the HH choice as 
we now show. An example of quantum cosmology with a Lorentzian region for small a was 
considered in the context of brane gravity in [25] . 

A generic boundary condition in region I that is not too far from a "stationary state" in 
the sense that it has comparable incident and reflected waves, will yield some finite ratio of 
B + to B-. Unless the coefficient of the increasing exponential B + is tuned specifically to 
vanish or to be much smaller than 5_ then the rising exponential will dominate the wave 
function inside the barrier and at the beginning of the Lorentzian region. Hence in practice, 
the HH boundary conditions can be taken, and will be a very good approximation to the 
generic situation. Of course, if for some reason (for instance, in analogy with a-decay) one 
would like to impose, as advocated by Vilenkin, that there is only an outgoing wave function 
in region III, then the coefficient B + would need to be very small. In our setup, this would 
seem to be a very special choice. 

The tunneling amplitude is of the form e ±iS? with 



$ = / da^J (a 2 — al_) {a\ — a 2 

J a- 



'(i-fti r^^fv (24) 

y k'sJ-n/WX o(t) 1 
In order to get the second equality we have made the substitution a 2 = a(r) 2 = (1 + 



1 — ^cos(2vAt))/(2A). A simple estimate of the integral which is quite sufficient for 
our purposes is obtained by assuming a triangular integrand ^JU(a) — v/\ 2 whose height 



is the maximal height of the potential barrier ^=wl — and whose width is a+ 



y^yl — y- The resulting estimate for $ is 

A more sophisticated analysis shows that the integral can be expressed in terms of complete 
elliptic integrals of the first and second kind. In the limit vj\ — > the exact result reduces, 
as it should, to the previous case (see eq. (jHJ)) and its value is 8/3, such that in this limit 
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$ = The exact expression is 



v/2 . 4za 1 1 f m — 1 r , 2 — m r - , 

$ = - — . — <^ 2 K\m\ + E[m] )■ , 26 

A v A ; / F 2m I m L J m L J < 



1 + 



where m = 2^1 - f /(l + yj 1 - f ). 

Let us now choose the HH + sign for the logarithm of the wave function for the reasons 
that were explained previously (we will discuss what happens when the other sign is chosen 
later). In this case, the tunneling amplitude is maximized at the maximum of $. Using the 
simple estimate of eq. (|25jh we find that when 

A = 8u (27) 

$ is maximized, and its value is 

®\max = (28) 

Using a more accurate numerical evaluation of the exact expression gives, 

A = 5.26z/, (29) 

and 

12 

<&\max = —■ (30) 

The important point here is that the ratio of the radiation energy density to the CC is 
some finite fixed numerical constant. This means that the initial radiation energy density 
determines the value of the CC, 

(3D 



7TC 4 Ml' 



It is clear from this formula that the result is sensitive to the initial temperature and that 
the precise estimation of the value of vj\ is not particularly important. This is essentially 
the same result as the one obtained in [11], except that they have effectively put the constant 
c = 1. Our derivation shows that setting c to unity is inconsistent, in that it would in effect 
give at the barrier a radiation energy density that is greater than the string energy density. 
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IV. THE TUNNELING PROBABILITY 

In order to discuss in a meaningful way the relative probabilities for tunneling into dif- 
ferent points in the landscape and the issue of whether or not inflation is favored one really 
needs, in addition to the CC, to introduce the set of moduli fields <fi and their potential 
V(4>). If the moduli potential is not steep, one could (following j23j) take over the results 
of sect. |H]with the substitution A — > 47r 2 cr 4 y(0). However, the string moduli potentials are 
steep except in a limited domain around their extrema. Let us ignore this for the moment 
and come back to this issue in section IVI1 

Let us consider the set of dynamical parameters in the potential that are determined by 
the fluxes, gauge groups, etc., and denote them collectively by (3. These vary from point to 
point in the landscape. Now the string to Planck mass ratio depends on the moduli ^f- (0). 
Using the above substitution, our previous maximization argument gives, 



^VX^^)^,^^^ . (32) 

It is clear from eq. ()32|) that the maximization of the tunneling amplitude puts a constraint 
on the parameters j3 and the values of the moduli fields to which the universe tunnels. 

Let us now observe what happens when the LV sign for the under the barrier wave function 
is chosen corresponding to the usual tunneling situation, with only outgoing waves in the final 
Lorentzian region. Then within the class of models that we consider: a closed universe with 
a positive CC, and radiation whose temperature does not exceed the Hagedorn temperature, 
the maximum of the tunneling amplitude becomes a minimum. Now the tunneling amplitude 
is maximized at the edge of parameter space when A = 4z/, exactly at the point that the 
barrier disappears! In this case eq. (J32j) should be replaced by 

i 



4 



4A4 U;^4,^g(^ . (33) 



Clearly, given that we are ignoring order one factors, the difference between the two cases 
is not that significant. The thermal boundary condition switches the physical consequences 
of the two wave functions. In the absence of radiation the LV wave function favors a larger 
CC whilst the HH wave function favors a zero CC. With radiation the HH wave function 
favors a larger CC than the LV. 

The maximization that we have performed was essentially with respect to A (or V) keeping 



v fixed. In 







probabilities for tunneling for different values of v are compared. However, 
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if we strictly follow their logic there appears to be a problem. From eq. (|3U|) we have, 



V2 y/2to?(*M. 



1 (34) 



I2u 12 n dof MY 

The number of light degrees freedom does not change that much - being around 10 2 — 10 3 
so the tunneling probability is essentially controlled by the ratio of the Planck scale to the 
string scale. In the heterotic string, for example, this is given by (See for example 2(| 
chapter 18). jfz = 9D ^ d - where a Y M = Qym/^ 2 i s the gauge field coupling strength a; the 



string scale. In type I theory on the other hand the ratio is given by (See for example 

m'i 

M) 



3)- 



g^^M. w here g is the string coupling jsjj]. Plugging these into eq. we get 



V2 9tt 2 c 4 8 

®max = -777 Heterotic, 

12 n dof a Y M 

V^QttV 4 

= 777 Type I. 

12 n dof ga Y M 

This seems to favor tunneling into very weakly coupled universes! If the value of M s /M p is 
not fixed then it is preferable to have it vanishingly small, and then the radiation energy 
density also vanishes and consequently also the CC. The final result in this case is very 
similar to the original HH result, the universe tunnels to the smallest possible value of the 
CC. We believe that to come to any reliable conclusion it is really necessary to explicitly 
consider the dependence on the moduli scalar fields, which we do in the next sections. 



V. ACTION AND WHEELER-DE WITT EQUATION FOR GRAVITY AND A 

SCALAR FIELD 

As we have explained, to determine the tunneling probability it is necessary to reexamine 
the wave function of the universe when in addition to the gravity sector, we also have scalar 
fields. For simplicity we will consider just one field with the action 

S * = ~ J d ' x V9 {^Td^A + V(<P)) . (35) 
For a closed universe the total mini-superspace action is 

S = i J dt | -ad 2 + a 2 a 3 <p 2 + a - An 2 a 4 a 3 V{<f>) } . (36) 
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Here we have absorbed the CC into the scalar potential V. The conjugate momenta are 



n a = —act 
n =aV 2 0, 



(37) 



and the Hamiltonian is 

H = -i^ + ^-l^-^ Wv ^- < 38 > 

The term n| determines the kinetic energy (KE) of the scalar field. 

The WDW equation is H^(a, 0) = 0, and can be obtain by the substitution II a — ► — ij^, 
and n^, — >■ — The result is then the following, 



a 2 n 2 - AfiJ + a 2 (a 2 - 47r 2 aVV(0)) 



*(a,0) = 0. 



(39) 



We have assumed a particular operator ordering but this is not particularly important for 
our arguments. 



A. A slowly rolling scalar field 

In the limit that the scalar field is moving very slowly and where 2a 4 V(<fi) can be treated 
as a constant, the WDW equation is separable \l/(a,0) = xi'P^i -)- Then, 



d 2 x(<f>) 

d<p 2 



-a 2 E X (<P) (40) 



a 2 ^^ - (a 2 - a 4 47rV\/(0) - ^(a) = 0. (41) 

The significance of is clear: it represents the average KE of the scalar field. The term 
E/a 2 in the WDW equation originates from a term E/a 6 in the energy density, which is 
indeed a scalar field KE term. 

Positive E corresponds to positive KE, and leads to an oscillating scalar field wave func- 
tion while negative E leads to an exponential scalar field wave function whose interpretation 
is unclear. The solutions of the scalar field WDW are simply linear combinations of "free 
particle" solutions: wave packets. For positive E, denoting k = \[E, 



X 



CO 

(0) = J dke lak *x(k). (42) 
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The choice E — leads to a constant wave function and a uniform probability for all values 
of the field. If one chooses E = 0, and the path integral is dominated by a single classical 
configuration, then this classical configuration is a constant scalar field, for consistency. 



If the potential is such that the field is moving significantly, we may use the following 
method to solve the WDW equation. We find the classical solution for the scale factor and 
scalar field, and express the scalar field as a function of the scale factor. This is possible 
provided that the scale factor is a monotonic function of time. We then use the parametric 
solution to find the kinetic and potential energy of the scalar field as a function of the scale 
factor as shown below. This results in a modified equation for the scale factor only. 

The energy density of the scalar field p^ = ^<p 2 + V(<j)) and its pressure is p^ = \4> 2 — V(4>). 
The conservation equation 



which is valid also for a closed universe, can be formally integrated if an additional relation 
between p and p is supplied, 



Let us consider as a simple instructive example the Euclidean scaling solution for which 
the total energy density of the field behaves as a fixed power of the scale factor: p^ = 



w^a) = p<f>(a) I p<$>{a) is constant. The constant C replaces the CC for the general case in 
which is not constant, and the solution is the same as the one with a fluid whose EOS 
parameter is w = p/p = m/3 — 1. In general, the scalar field will not have a constant EOS 
parameter. 

In the presence of space curvature it is difficult to find potentials that yield a constant 
EOS. For the Lorentzian spatially flat case it is well known that exponential potentials lead 
to scaling solutions. A potential of the form V((p) = Ae a ^ leads to solutions that have 
4> 2 ~ V(4>) ~ l/t 2 and w = — 1 + 2a 2 /3. For a < 1 the potential is flat enough and the 
expansion is accelerated. The maximal EOS parameter for a scalar field with a positive 
potential is w = 1 in the limit that the evolution is KE dominated. This limit is reached 



B. A rolling scalar field 




(43) 




(44) 



|0 2 + V((j)) ~ C/a m , corresponding to the case that the equation of state (EOS) parameter 
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when a approaches the critical value v3, and A vanishes. For a scaling solution to exist 
for larger values of a, the prefactor A needs to be negative, and in this case values of the 
EOS parameter that are larger than unity are allowed. A positive potential that is steeper 
than the critical steepness leads very quickly to a KE dominated evolution, as discussed for 
example, in [28]. This will be relevant to our discussion in sect. I VI II 

The Euclidean equations are in fact the same as the Lorentzian equations for the spatially 
flat case, and therefore the Euclidean scaling solutions lead to the same p^{a) as do the 
Lorentzian ones. 

In the WDW equation (j3HJ), the terms that correspond to p^ are — 2 J a2 II| — 2a 4 a 4 V r (0), 
so on a scaling solution we obtain the following effective WDW equation for the scale factor 
only, 

(a) = 0. (45) 



r\2 

+ a 2 -Ca 4 ~ m 



da 2 

One can immediately realize that for the potential barrier to exist and be finite, one needs 
that m < 2, so that the scaling behavior corresponds to a power law inflation with an 
equation EOS p/p that is more negative than that of spatial curvature p/p < —1/3. In 
general, on a solution p^ can be replaced by some function of a. The potential is then 
U(a) = a 2 — Pcj>{a)- Clearly p^ needs to grow at large a faster than a 2 , and therefore needs 
to have an EOS more negative than that of curvature. 

Recall that the EOS of space curvature is —1/3 which is just the borderline between 
decelerated and accelerated expansion. Thus, our conclusion is that the universe tunnels to 
a state of accelerated expansion. The scalar field energy density needs to dominate over the 
space curvature for a while, in order to get the potential down to zero, where the Lorentzian 
era begins. This conclusion will be very important for the following discussion about the 
landscape. 

An additional conclusion is that if the dominant power in p^ has a less negative EOS, 
the larger the range of a that is required for the "dark energy" to become dominant over 
the space curvature. This is depicted in Fig. ^ In addition, for lower dominant power in 
P0, the faster the field <j) is moving. As the evolution becomes less and less similar to that 
in the presence of a CC, the scale factor is moving more when under the barrier and the 
scalar field needs to move more. This requires that the potential of the scalar field be flat 
over a larger region in field space so the accelerated expansion could be supported. Another 
general comment is that one can compare the tunneling amplitude for a CC to that in the 
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U(a) 




a 



Figure 1: The potential of the WDW equation for a closed universe with thermal radiation. The 
solid line corresponds to the potential for a universe with a CC, and the dashed line corresponds 
to a universe with dark energy with a softer EOS. Regions I, III are classically allowed regions 
where the wave function is oscillatory, while region II is a classically forbidden region where the 
wave function is exponential. 

case of "softer" EOS (less negative than —1). Assuming that the EOS is a constant, leading 
to a power law dependence p^ ~ C/a^ 2 ~ m \m > 2, we can see that the tunneling amplitude 
that comes form the positive exponential mode (in the case that the magnitudes of the 
potentials are similar) is much larger for the softer EOS, since the barrier is much larger. 
This is depicted in Fig. ^ 

It should be stressed here that the considerations of this section are independent of the 
specific choice of the boundary conditions on the wave function. We have argued that 
generic initial conditions lead to the dominance of the rising mode under the barrier. If 
however one imposes final conditions such that only the outgoing wave is allowed, then the 
usual tunneling picture emerges. Obviously, in this case too the barrier needs to end and 
therefore, as argued above, the tunneling needs to take place to an accelerating universe. 
By altering the topology of the final state as in [22] (or even just the geometry by taking 
flat or negatively curved universes) one could eliminate the barrier altogether and in this 
case our arguments would be irrelevant. 



18 



VI. THE TUNNELING PROBABILITY WITH MODULI 

The discussion and results of sections IHII and IIVI about the tunneling probability and 
its extremization did not take into account the fact that the ratio of the string scale to 
the Planck scale M s /M p is dynamical and moduli dependent, and that the CC should be 
replaced by the moduli potential. We reevaluate them in the light of the discussion of section 
IVl We make several assumptions about the form of the moduli potentials that allow us to 
obtain more definite results. 

In the case that the scalar field is a modulus, there are two possibilities. First, that 
the true CC is large and then the potential is dominated by the CC. This case is not very 
interesting because it does not lead to a universe that is similar to ours. The more interesting 
possibility is that the true CC is substantially less than the string scale. In the latter case 
we know from general arguments about Peccei-Quinn symmetries and how they break, that 
the potential is a sum of steep functions. Further, we know that at generic points that are 
not extrema of the potential and that are in the outer region of moduli space, the potential 
is dominated by a single steep function (see, for example, a related discussion in j^). It 
follows that the only flat regions in the potential of the moduli fields where V'/V <C 1 occur 
in the vicinity of an extremum. 

We argue that consequently, when the CC is small the universe cannot tunnel to a generic 
point in the outer region of moduli space. We have seen that the tunneling end point needs 
to be where the potential energy dominates the energy balance in the universe, and in 
particular dominates the field's KE. However, we have just argued that the potential for a 
generic point in the outer region of moduli space is steep. On the other hand, we know that 
the cosmological solution of a scalar field on a steep potential leads to the dominance of KE 
over the potential energy. Our conclusion follows. Additionally, we have just argued that 
the only flat regions in the potential of the moduli fields are in the vicinity of an extremum. 
Our tentative conclusion is therefore that the universe tunnels to a region that is not far 
from an extremum of the moduli potential. We proceed to examine this conclusion in a 
more concrete setup. 

A possible loophole in our argument could exist if some hitherto unknown Euclidean 
transient solutions that support accelerated expansion for a short period of time could 
be found. We expect such solutions, if they exists at all, to require some special initial 
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conditions. We cannot discuss their possible existence in the general setup without adding 
more specific information on the moduli potentials. In sect. lVIll we show that such solutions, 
even if they do exist, do not modify our conclusions. 

The moduli potential can be put in the form of an N = 1 SUGRA potential, 



V = e K 

with 



DiWDiWK* -3\W\ 



(46) 



K = -]n(2S R )-3\n(2T R ) + ... (47) 

Here Tr is the real part of T and Sr is the real part of S and the ellipses represent the con- 
tribution of the other (Kaehler and complex structure) moduli. The field S is the (complex) 
dilaton axion field whose expectation value determines the coupling and T is the so-called 
volume modulus whose expectation value determines the size of the internal manifold. To- 
gether, they determine the string to Planck mass ratio Mf/Mp = S^T^ 3 . Let us denote 
the set of dimensionless moduli (S,T, ...) by <ft l . The potential is a function of these moduli 
as well as of the flux parameters and Casimirs of the gauge groups which we collectively 
denoted by (3 as in the discussion above eq. (J52j) . Thus we may express the potential as 



V = V(<F;I3). (48) 



The WDW equation in the presence of the radiation term becomes, 

-a* + K*-^ + a 2 (a 2 - a 4 V(0; (3) - K) 

\daj dfrdft V K ' ' 



#(a,0) = O. (49) 



Here we have substituted A — >• V as in sect. |V| We can reason as in sect. IV Bl that tunneling 
to an accelerated expansion phase is favored since otherwise the barrier does not end. Since 
an accelerated expansion phase requires that the potential is flat enough, it reasonable 
to expand to leading order in an expansion in V'/V, and as argued previously it is also 
reasonable to assume that the moduli are near a critical point in the potential. Thus the 
arguments of sect. II I II that lead to eq. (jS2} can be used. The new probability maximization 
conditions are 

V(fr/3) ~ 5z/(0), (51) 
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with 

n d0 f Mj _ ^L( T) -s {Sr) -i_ (52) 



v 



9tt 2 c 4 M 4 9tt 2 c 4 



Equations (f5T)|) (one for each modulus) determine, for any set of parameters f3, a region in 
moduli space for which the slow roll conditions are satisfied. Of course, there may be some 
values of f3 for which there is no solution. Equation (JoTj) is then a further constraint that 
restricts the parameter values only to those that satisfy it. 

Equations ()50)) - ([52j) . however, leave the system under-determined. The only exception is 
when there is only one parameter in the set (3. The logarithm of the WKB wave function 
at the point of emergence from the barrier into the Lorentzian spacetime depends on the 



values of the moduli at this point {S L ,T L , z L } = {<fif} and on a 



+ ■ 



(53) 



Hence, it takes the following form in general 

^ = ^(S L ,T L ,z L ,a + ,P). (54) 

The wave function should be extremized with respect to all the moduli as well as the pa- 
rameters p. Since $ depends on (3 only through its dependence on the potential, we have 

d<& <9$ dV 



9(3 dV\tfdp\tf 
<9<3> <9$ da+ <9$ 



0, (55) 
= 0. (56) 



d<p l L da+ltfdtf d(j) l L \ a+ 

In principle, these equations determine the tunneling end-point in the moduli space that 
serves as the initial values for the classical evolution of all the moduli. They also determine 
the discrete parameters, the fluxes, gauge group parameters, etc.. These equations therefore 
determine the particular flux configuration to which the tunneling occurs. Of course, in 
general there may be more than one solution to these conditions, so their solution may be 
a multiverse rather than a universe. 

In practice we used the particular expression for $ that was obtained in eq. (|26|) for the 
case of a constant potential (CC) by arguing that the tunneling took place to a flat point 
on the potential. Thus we argued that even with the moduli, the equations of sect. |^ 
could be used with the appropriate replacements for A and v. Equation (}5T|) then follows 
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from imposing 8<&/dV = which implies that the whole set of equations ([55)1 . one for each 
parameter in the set (3, are satisfied. Clearly this is a sufficient but not a necessary condition 
for extremization with respect to (3. It would be interesting to explore the nature of the 
general solutions to the extremization equations. 

In the next section we will investigate the nature of the restrictions that we have found 
for the specific potentials for moduli that have been suggested recently by KKLT. 



VII. APPLICATION OF THE THERMAL BOUNDARY CONDITION TO THE 

KKLT MODEL 

We would like to determine more accurately the point on the moduli potential to which 
the universe tunnels. For this we need to input some additional information on the properties 
of the moduli potential. We will use here the KKLT model In this model, the potential 
has two contributions Vsugra and Vjj. The first takes the form of an iV = 1 SUGRA 
potential (using the same units as before) with a Kaehler potential given by eq. ()47j) and a 
superpotential 

W = A + BS + Ce~ aT . (57) 

The first two terms in eq. (|57j) come from the fluxes Q|, A and B are functions of the 
complex structure moduli, and there is only one Kaehler modulus T. The third term in 
this expression can arise from gaugino condensation in a gauge group living on a stack of 
7-branes wrapping a four cycle on the compact manifold. KKLT assume that it is possible 
to ignore the third term and integrate out S and the complex structure moduli, assuming 
a flux configuration which makes their masses heavy. While this is not strictly correct (see 
[30! ]) the corrections are not important for the current discussion so we will ignore them. 
Then S and the complex structure moduli are constants, and the effective superpotential is 
of the form W = Wq + Ce~ aT with a Kaehler potential K = —3 ln(T + T) giving 

V SUGRA = ^fT~ \W* + ^ T * a + ^Ce~ aT -j . (58) 

The minimum of Vsugra is a t DW = dxW + dxKW = 0, and for consistency with the 
assumption that the volume is large, and therefore that ten dimensional supergravity is 
valid, and for consistency with the expansion in non-perturbative terms, one needs Tr ^> 
1, aT > 1, so that we need to have W < 1. This can be achieved by fine tuning the fluxes. 
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The second contribution to the KKLT potential Vjy comes from anti--D 3 branes and 
breaks supersymmetry explicitly from the 4D perspective. It is proportional to ^ in a 
naive calculation of the anti-Dbrane tension, however, a term proportional to has also 
been proposed in the literature. The total potential has a shallow positive minimum at a 
largish value of Tr, and a small barrier separates it from a positive tail that goes to zero as 
~ 1/21. 

Now let us check whether the end point of the tunneling can be on the asymptotic 
tail using the arguments of section IV Bl The field Tr is not canonically normalized. The 
canonically normalized field is x where Tr = . So the asymptotic dependence of the 
full potential on x is V ~ e~ 3 ^ x . Recall the discussion in sect. IVI where we concluded 
that the (Euclidean) cosmological scaling solution of a canonically normalized scalar field <fi 
with an exponential potential V = Ae a ^ gives a power law dependence for the scale factor 
a(t) ~ t Pa , with p a = 1/a 2 . If p a > 1, a < 1, the expansion is accelerated. Recall also that 
for \a\ > y/3 the prefactor A needs to be negative for a scaling solution to exist. In the case 
that the potential goes like 1/Tr then a = a/6 or in the case that the potential goes like 
1/2r then a = -^/|, both significantly above I. The prefactor d is positive. Hence we can 
conclude that these potentials do no lead to accelerated expansion. 

The argument of section IV Bl indicated that for the tunneling barrier to be finite, one 
needed an accelerating scale factor. Our conclusion is therefore that the tunneling end point 
cannot be on the asymptotic tail of the potential in the region Tr 3> 1 where Vjy dominates. 
A possible loophole in the argument is that perhaps it is possible to find a transient solution 
that includes a brief period of accelerated expansion. However, also in this case the potential 
energy has to dominate the energy balance. Then, ignoring numerical factors of order unity 
and since Sr is fixed at a number of 0(1), eq. lJoTj) leads to the condition 



9vr 2 c 4 



d. (59) 



As we discuss below, this relationship cannot be satisfied since d has to be very small, and 
the l.h.s. is not particularly small. Hence it is not possible to satisfy the extremization 
conditions even in this case. 

If the anti-Dbrane term is not included then the tail of the potential is steeper, and if a 
potential of the form d/T^ is added then eq. (Jol^ is replaced by = dT R which is even 
harder satisfy. Our conclusion is therefore valid for these additional cases. 
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Let us now check whether the tunneling end point can be near the shallow minimum, or 
near the maximum of the barrier separating the minimum from the asymptotic region. The 
dimensionless radiation v is given by eq. ((H2} and the dimensionless potential by 

V = Vsugra + ^- (60) 

R 

At the two extrema, the two terms are comparable. In the event that the tunneling occurs 
to an extremum, which we have argued is a reasonable approximation, the wave function is 
extremized at points where v and V are related as in eq. (|51|) . In our case Sr fixed at a 
number of 0(1). 

Equation (joTj) gives 

5^ = aAT^ ^ + ( l TRa + 1)Ce _ aT ^ + ^ (61) 

where we have ignored a factor of 0(1). The factor is about O(10 _1 ), and is not 

expected to vary much with different flux choices. On the other hand, for models with one 
condensate (i.e. one non-perturbative term as in the original KKLT example) and a small 
CC, both extrema are at quite large values of clTr. For large values of clTr the r.h.s. of 
eq.(jnU is much smaller than 1/10. We conclude that eq. lJoT]) cannot be satisfied at either 
extremum. For example, let us consider the choice of parameters in KKLT, where a = .1, 
A — 1, Tr ~ 100 — 150 and d — 1 x 10~ 9 . Although the normalization of the potential in 
KKLT is somewhat different than ours, it is clear that the r.h.s of eq. (jfilj) is of order 10 -9 
for the maximum, and much smaller for the minimum where the CC was tuned to be very 
small. 

Our conclusion is that the quantum cosmology arguments would prefer a modified KKLT 
model, with additional exponential terms in the superpotential to allow extrema that are 
closer to the central region of moduli space where eq. (joTj) has a chance to be obeyed. Such 



models can be constructed 
(see, for example 



3V having several non-perturbative terms in the superpotential 
I]). Clearly, it is easier to satisfy eq. (JSTj) at a maximum, a 
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saddle point, or a metastable minimum rather than at a global minimum. 



VIII. CONCLUSIONS 



We have seen that the probability distribution obtained from the wave function of the 
universe can provide an interesting and restrictive dynamical selection principle on the land- 
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scape of string solutions without reference to the anthropic principle. 

We have proposed that the universe emerged from the string era in a thermally excited 
state above the HH vacuum, and determined, within a mini-superspace model, the proba- 
bility of tunneling to different points on the landscape. We have clarified the significance of 
including a radiation term for the HH wave function of the universe and have shown that the 
radiation term switches the roles of the Hartle-Hawking and Linde-Vilenkin wave functions. 
We have found that the potential energy of the tunneling end point from which the universe 
emerges and begins its classical evolution, is determined by the primordial temperature, and 
that this starting point can be followed by some interesting dynamics. 

We have found that a more accurate treatment, even within the mini-superspace approx- 
imation, requires the inclusion of the moduli fields, and we have included them. By taking 
into account some generic properties of the moduli potential we then argue that the tunnel- 
ing to the tail of the moduli potentials is disfavored, that the most likely emergence point is 
near an extremum, and that this extremum is not likely to be in the outer region of moduli 
space where the compact volume is very large and the string coupling very weak. Combined, 
these considerations select a class of values of the flux parameters etc., that characterize a 
universe, or a multiverse, in the landscape. 

We explicitly demonstrated the applicability of our arguments for the KKLT model of 
moduli stabilization. We have determined that for the KKLT model the tunneling to the 
tail of the potential or to the vicinity of the barrier that separates the minimum from the 
asymptotic region is disfavored. Our quantum cosmology arguments favor tunneling to an 
extremum of the potential that is close to the central region of moduli space as might be 
obtained from generalizations of the original KKLT model. 

Finally, we might consider the relevance of the counting program of Douglas and collab- 
orators a El to our arguments. We have calculated a quantity that is analogous to the 
square of the tunneling amplitude in the calculation of a decay rate of an unstable parti- 
cle. To obtain the total tunneling rate one needs also the density of final states. Perhaps 
the counting program can supply the latter, so that a complete calculation of the relative 
probabilities of finding one or another universe, or a certain subset of the multiverse, is 
obtained by taking the product of the two factors. The number density of solutions does not 
seem to influence much the preferred value of the CC, because the tunneling amplitude is 
sharply peaked as a function of the CC, while the number density of solutions is uniform as 
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a function of the CC. Perhaps it is more relevant to the preferred values of other dynamical 
parameters. 
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Appendix 



Here we explain the relation of our work to the calculations of Sarangi and Tye 11 1. 

Sarangi and Tye (ST) have argued that the wave function needs to be modified due to 
decoherence effects. They argue that the fluctuations of the metric and of other light fields 
should be integrated out and traced over. Their calculation is rather involved, but the final 
result can be justified and related to that in sect. IHII bv the following simple argument. 

ST are essentially computing the thermal partition function in a Friedman-Robertson- 
Walker (FRW) background at some undefined temperature. It is well known that this 



ree energy) which is just the energy 
I. It is also well known that in the 



calculation gives a contribution to the effective action 
density of radiation i.e. a term proportional to T 4 [i; 
FRW background the temperature scales as the inverse of the scale factor a -1 . So the only 
question is: what is the temperature? 

Let us work with the action Q. The effect of the unobserved fluctuations can be repre- 
sented by the entropy S within a horizon volume. An order of magnitude estimate of the 
upper bound on this is (a related though not identical calculation is performed by ST) 



d n 



H 



h J \VAl s 

where we have assumed as in |ll( that the infrared cutoff is the horizon size and the ul- 
traviolet one is the string scale l s . Also we have used the Friedman equation to estimate 
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H ~ 1/ \/A. Comparing to the entropy of a thermal state 5 1 ~ a 3 a 3 T 3 , we expect S" and S to 
agree to within a time independent factor c 3 . The factor c may be numerically large but is 
not expected to be parametrically large. More importantly it is time independent. Equating 
them ^—^j-^j = c 3 a 3 a 3 T 3 we find that the effective temperature is T = l/(ca^/\l s ) where 
A is the dimensionless CC introduced after eq. (0). The energy density associated with this 
radiation is then given by eq. (fT3|) . ST effectively have c = 1. However this is not consis- 
tent with the effective field theory mini-superspace starting point, since as can be seen from 
eq. (H3), at a — a± ~ 1/ (A^), this would give a radiation energy density that is greater than 
string scale. So for consistency one should have c 4 > n do f as discussed in sect. IHII 
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Note that in our universe oym — 1/25 so given that M p ~ 10 19 GeU we get a string scale of 
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about 10 GeV if the heterotic string is the correct theory whilst a similar result is obtained 
if g ~ 0(1) in the type I case too. This is the value we quoted earlier in the discussion on 
inflation. 

[35] ST effectively work with a thermal circle whose radius is called T. This should in fact be 
identified with the inverse temperature - so in our notation their T should be replaced by our 



